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^ We study the problem of designing electrode structures that allow pairs of ions to be 

. brought together and separated rapidly in an array of linear Paul traps. We show that 

' it is desirable for the electrode structure to produce a d.c. octupole moment with an 

' a.c. radial quadrupole. For the case where electrical breakdown limits the voltages 

^ that can be applied, we show that the octupole is more demanding than the quadrupole 

, when the characteristic distance scale of the structure is larger than 1 to 10 microns (for 

' typical materials). We present a variety of approaches and optimizations of structures 

' ■ ^ ' consisting of one to three layers of electrodes. The three-layer structures allow the fastest 

^ operation at given distance p from the trap centres to the nearest electrode surface, but 

, when the total thickness ui of the structure is constrained, leading to «) < p, then two- 

' layer structures may be preferable. 
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^ . ■ 1 Introduction 

^ ' The application of ion traps to quantum computing El El has lead to interest in the 

. ^ , construction of systems of multiple ion traps situated close together. Experiments have shown 

' that quantum logic gates with good fidelity can be produced between a single pair of ions in 

i a single trap, and such methods can be extended to a few ions in a single trap El El EUHl . 

However, to scale this approach up to the manipulation of many qubits, it is not feasible to 
try to manipulate large numbers of qubits in a single trap. Instead, it is desirable to have 
an array of traps, with the possibility to move quantum information between ions in separate 
traps. A method to move the quantum information which has advantages of relative simplicity 
and robustness is simply to move the ions themselves [51 HI ITU]. To achieve this we need an 
array of ion traps, with the possibility to move ions between traps. This concept has been 
demonstrated in an impressive set of experiments I12L I13| which demonstrate its promise. 
The present work addresses the issue of how to design a basic element of such an array: a 
region in which a pair of ions can be brought together or separated rapidly, by adjusting d.c. 
voltages on electrodes. This includes the case of separating into two wells a pair of ions which 
are initially in a single harmonic well. 
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There are further constraints which make this a non-trivial problem: 

• The electrode structure should allow a route for ions to be moved into and out of the 
region of the close pair of trapping centres. 

• The surfaces of the electrodes should be kept as far as possible from the trap centres 
where the ions sit. This is in order both to minimise heating of the motional state of 
the trapped ions (see below) and to reduce the impact of surface irregularities on the 
electric potential fimction experienced by the ions. 

• The trap confinement should be tight, i.e. the vibrational frequency v of the ions must 
exceed some desired minimum value. This is in order to maximise the speed at which 
ions can be moved by displacing the trap centres. 

• The electric fields at the electrode surfaces must not exceed the electrical breakdown 
limitations of the materials. 

• The electrode configuration must be capable of being fabricated accurately at the re- 
quired distance scale. This may require microfabrication techniques which have their 
own intrinsic limitations. 

For diagnostic purposes it may also be useful for the electrode structure to be open in 
order to allow optical access. 

The need to keep the electrode surfaces far away from the ions is largely owing to the 
observation of an anomalous heating rate in ion traps, which becomes a severe problem with 
traps of small distance scale. The heating process is not fully understood but appears to 
be associated with impurities deposited on the electrodes, and studies suggest the heating 
rate scales as where p is the distance from trap centre to the nearest electrode surface 
|14[I11[[TH] . Electrodes further from the trap centre can also be bigger, which usually implies 
they are easier to make. 

Typical distance scales are 1-10 microns for the separation of the trapped ions, 1-100 
microns for the distance to the electrodes [HI ITTl [HI ITH] . 

We begin the discussion by considering in section |21 the general problem of a pair of 1- 
dimensional potential wells which can be moved together or pulled apart while maintaining 
the tightness, i.e. the normal mode frequencies of a pair of ions trapped in the wells. We 
argue that the best electrode structure is one which produces an electric octupole potential, 
which is then 'tweaked' by the addition of quadrupole terms. In section|21we introduce radial 
confinement by an a.c. quadrupole, and discuss in general terms the expected scaling of the 
trap parameters with the size of the electrode structure. We show that for structures whose 
dimensions are greater than 1 to 10 microns, the most demanding factor in the design is to 
obtain a large octupole moment, assuming the applied voltages are limited by electrical break- 
down of the materials. We then proceed in section^ to consider the problem of obtaining an 
electric octupole with a minimum number of electrodes. In section [S] we describe numerical 
calculations of the electric potential for several types of electrode structure, obtaining the 
values of two dimensionless factors that give the strength of the three-dimensional d.c. oc- 
tupole and two-dimensional a.c. quadrupole moment at given breakdown field and trap scale. 
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Section El then briefly analyses the influence of manufacturing imprecision which breaks the 
symmetries of an electrode structure and introduces stray electric fields. 

Section 13 discusses the main features of the numerical results. For example, we quantify 
the relative merits of a 3-layer over a 2-layer or single-layer (planar) design. We estimate 
the impact of requiring the complete structure to be thin (to aid microfabrication). We sum- 
marize the optimizations obtained by adjusting the positions and relative sizes of electrodes. 
Section IHl concludes. The appendices give a discussion of axial micromotion and some general 
information on charge distributions that produce electric octupoles. 

2 Generic study of V{x, y, z) for two traps close together 

We consider first of all some general properties of the electric potential for two Paul traps 
located at {x,y,z) = (0,0, ±s). We assume the vibrational frequencies , cjj, , cj^ of a single 
trapped ion are the same for the two traps. The precise relative sizes of lOx, ujy and ujz are 
not crucial, but it will be useful if LOx^LOy ^ oJz- The electrode configurations to be discussed 
in later sections of the paper will produce the confinement along z primarily by d.c. voltages, 
and that along x and y primarily by an oscillating (r.f.) quadrupole field. However, we do not 
need to assume this for the general discussion in this section (for clarification of this point, 
see the appendix). 

Let the distance from the origin to the nearest electrode surface be p. Since we aim to 
achieve close traps and far electrodes, we assume s p. We can then usefully analyse the 
electric potential V near the origin by a Taylor or a multipole expansion. We further assume 
that the electrode construction is close to symmetric so that odd order terms in the expansion 
almost vanish. First let us consider V along the z axis: 

V{0,0, z) Vo - Eqz + az^ + (3z^ (1) 

where we have dropped the cubic term since we assume that it is negligible. The signs of the 
coefficients Eo,a,(3 have been chosen so that Eq is the electric field at the origin (assuming 
dV/dx — dV/dy — there), and to create a double well potential we need a < and /3 > 
(see figure 1). 
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Fig. 1. Plot of Viz) (equation Q) for Vq = 0, Eq = 0, a = -1, 13 = 1. 

First suppose there is no linear term, Eq = 0. By solving dV/dz = we find for a < 
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the potential wells are centred at 



s = 




(2) 



Evaluating d'^V/dz'^ at either well centre gives d'^V/dz'^ — Aa. If the oscillating voltages are 
so arranged that the z axis is an r.f. null, i.e. V{0, 0, z) is constant in time or only has a 
relatively small oscillation, then we may obtain for an ion trapped in either well directly 
from d^V/dz"^, giving 



where m,q are the mass and charge of the ion. For example, at lo^jl-K = 1 MHz, \a\ ~ 4 
MV/m^ for the charge to mass ratio of the ^■^Ca+ ion. The well depth V^(0) — V{s) = 
e.g. 0.2 mV (2 Kelvin) at s = 10 ^m. 

In summary, the quadratic term a in the Taylor expansion sets the trap tightness (leading 
to LOz), and the ratio of the quadratic to quartic terms sets the trap separation. 

We will arrange for the linear term Eg to be small by using electrodes symmetric under 
reflections in the x—y plane, and then nulling any remaining field if necessary by compensation 
voltages. Let us calculate next the value of Eq which "tips the potential over" just enough 
so that one of the traps cannot confine ions, i.e. d^V/dz^ goes to zero at a trap centre. This 
happens when 



For example, for '^^Ca"'" with cj/27r = 1 MHz and s = 10 /im, \a\s ~ 40 V/m. 

2. 1 Cancellation of quadratic term 

Consider the Taylor expansion for a set of electrodes of distance scale a. In the absence 
of special choices or symmetries which make one or more terms disappear, we would obtain 



(for example, consider the electric potential due to a point charge, or a line charge, at distance 
a from the point or line). Therefore, without a special design, the order of magnitude of s 
is expected to be s ^ a. To obtain s ^ a we require a/P <C a^, and to obtain this without 
reducing lUz we must increase /3. This can be done by constructing an electrode configuration 
in which the quadratic term a is suppressed compared to /3, and then increasing the voltages 
on all electrodes together. Assuming the increased voltages are attainable (and we will find 
they are), the available (3 at given a is limited by electrical breakdown associated with high 
electric fields at the electrode surfaces. 

Assume that the electrodes have reflection symmetry in x-y, x-z and y-z planes, so that the 
electric field vanishes at the origin. Then odd-order multipole moments of V{x, y, z) vanish, 
and so do mixed 2nd derivatives such as d^V j dxdy at the origin. Under this assumption. 




(4) 



\E\^ , |a| ~ — , \(3\ 

a a 




(5) 
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we will next show that the required electrode configuration is one which produces, at lowest 
order in its multipole expansion around the origin, an octupole moment. 

An electrode set that can produce a/ (3 <C is close to the condition a ^ 0. Since we 
require the possibility to adjust a and P (e.g. by adjusting voltages on the electrodes) the same 
electrode set can also attain a = (in practice, it may or may not actually be used in that 
condition). We therefore have to consider electric potential configurations having d'^V/dz^ = 
and d'^V/dz'^ 7^ 0. From Laplace's equation, we then have d^V/dy^ = —d^V/dx^ , thus a 2D 
quadrupole moment in the x — y plane unless these derivatives are also zero. Unless we take 
special measures to make \d'^V/dx'^\ small, we have \d'^V/dx'^\ ~ /3a^ 3> muuHq when s <C a. 
This implies the field is strongly expelling along some direction in the x — y plane. To permit 
the ions to remain trapped, either the whole field must be made to oscillate, or an oscillating 
quadrupole field must be added whose strength is sufiicient to overcome this expulsion. If 
the whole field oscillates then in the pseudo-potential model, the effective values of both a 
and (3 are reduced, which is counter-productive. If instead a separate oscillating quadrupole 
is added in the x — y plane (with small effect in the z-direction) then its strength must be 
large, and this leads to electrical breakdown which will constrain the available range of (3. 

We conclude that the d.c. electrode structure must be designed to produce small \d'^V/dx'^\ 
(as well as small \d'^V / dz'^\) and therefore we require an electrode structure at or close to an 
electric hexapole or octupole configuration. A quantitative statement of this condition is given 
in sections 13 (eq. (|21|) ') and 13.21 The octupole has the advantage that refiection symmetries 
which cancel odd-order multipole moments can be used to avoid producing an unwanted 
electric field (i.e. a contribution to the electric potential varying linearly with distance at the 
origin), and therefore we will concentrate on that case. 

Confinement of the ions in three dimensions is completed by adding to the d.c. octupole 
field an oscillating 2D quadrupole field whose size is controlled separately. Hence, after 
reinstating a non-zero value for a, we obtain the general form 



where a, (3 and Qac are time-independent, Vi{x, y, z) is an octupole potential with 14(0, 0, 2) = 
z^, and is the frequency of the a.c. quadrupole providing confinement in the radial direction. 
Many ion trap configurations in common use are described by ©, but we have in mind a 
case where f3/\a\ is as large as possible. We will see later that none of the electrode sets to be 
discussed realise © exactly, because the time-dependent part depends slightly on z as well 
as X and y, and because higher-order terms appear, but it is useful to make clear what we are 
aiming at. 

2.2 Coulomb repulsion and normal mode frequencies 

If there is one ion in each trap, then the distance between ions is not 2s owing to their 
mutual Coulomb repulsion. The equilibrium positions are z = ±d/2 where the separation d 
is a real positive solution of 




-l-QacCos(m)(a;2 -2/2) 



(6) 



(3d^ -f 2ad^ = 

27reo 



(7) 
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In order to bring out the general behaviour of solutions to this equation, it useful to express 
it in the form 

where 

47reo|a|(2s)3- ^''^ 

Therefore the solutions are determined by the sign of a and a single parameter e, which 
compares the Coulomb repulsion force to the trapping force for ions situated at ±s in a 
harmonic well. When the /? term in Q dominates, e is large, and d ~ ((7/27reo/3)^^^- When a 
is large and positive, i.e. two ions in the same harmonic well, d ~ (q/Aneoay^^ . When a < 
and e <C 1, i.e. a pair of well-separated traps, then d ~ s(2 + e). 

The system of two ions has two normal modes for motion along the z direction: the centre 
of mass mode of frequency uJi and the breathing (also called stretch) mode of frequency UJ2- 
These frequencies are given by 

Lol = {2a + Sf3d'^)q/m, (10) 
L^l = w?(l + e), (11) 



where 



TTeomuJic 



(12) 



(to was defined previously, it is the mass of a single ion). 

When the ions are far apart, e — > e ^ 1 and the two mode frequencies are both equal to 
the oscillation frequency of a single ion in either trap. When a = 0, e = 2/3 and when a 3> 1 
(i.e. a single harmonic well), e — > 2. 

The reason to bring the ions close together is in order to perform a two-ion quantum logic 
gate in which their Coulomb repulsion provides the interaction energy between the qubits. 
The gate is fast when the mode frequencies are high and well-separated jl7l I18| . Therefore 
we would like a high value of both lui and e. Suppose the ions are initially in traps far apart, 
so that e ^ 1. Then the vibrational frequencies are almost equal and are given by the value 
of a through cq. 0. To bring the ions together, the voltages on the electrodes may be 
changed so that /3 gets larger while a stays the same size. Eventually electrical breakdown 
prevents further increase of If at this point the values of lui and e are high enough to 
give an acceptable gate speed, and if it is possible to apply the logic-gate laser pulses to ions 
in the separation zone, then the gate is performed and the ions are subsequently separated. 
Otherwise, to bring the ions into a common trap, |a| must be lowered and one accepts an 
unavoidable reduction in mode frequencies while the potential barrier between the twin traps is 
lowered, see figure|21 After a passes through zero and becomes positive, the mode frequencies 
increase again and e rises towards its maximum value of 2. The ions can then be moved as 
a pair to a convenient (e.g. less noisy) region of the complete array, for the logic-gate laser 
pulses to be applied. 
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Fig. 2. Plot of the centre-of-mass frequency uii (solid curve) and stretch mode frequency (dotted 
curve) vs. a for three different values of f3. The curves shown are for /3 = 100, 10, 1 in arbitrary 
units, the higher value of /3 giving the higher normal mode frequencies. The figure shows the 
generic behaviour, the units are arbitrary. 

At the condition a — 0, the fuU analysis is tractable and one finds 



and e = 2/3 as mentioned above. 

In the rest of this paper we will consider methods to produce this condition, i.e. an 
octupole potential, with an added oscillating quadrupole potential in the x — y plane, as in 
(0. In practice the electrode structures to be considered depart from @ primarily in one 
respect, namely that the oscillating term includes a z-dependence, i.e. 



where we introduce the parameters ax^oiy^ to characterise the 3-dimensional a.c. quadrupole. 
If 7^ ay then necessarily 0, and in this case there is axial micromotion. It is shown in 
the appendix that this leads to a Mathieu equation, which can be solved in the standard way 
as a combination of slow secular motion and fast micromotion [201 121| . The secular motion 
can be modelled as motion in a pseudo-potential which takes into account the influence of 
both the d.c. and the a.c. (ponderomotive) terms. The above analysis then remains correct 
(within the pseudopotential approximation) as long as we replace a by a', where a' is the 
coefficient of the term in the total effective potential. For the structures discussed in this 
paper, it is found that the electrode voltages required to produce a' — are very close to 
those required for a = 0, hence changing from one condition to the other has little effect on 
(3. At a' = the secular frequency is completely determined by the value of /3 through H13(l . 




(13) 




(14) 



3 Radial r.f. confinement and scaling 

Consider the potential A pair of trapped ions is confined in the radial (x and y) 
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direction by the oscillating quadrupole field. To first approximation we may ignore the effect 
of the octupole potential on the radial motion, then this motion is described by a Mathieu 
equation. In the pseudo-potential model, for a = 0, the secular frequency of the centre-of-mass 
vibrational motion in the radial direction is given by |19l 1201 12] 

uJr ^ qrn/2V2, (15) 

where 

is the Mathieu g-parameter. For stable motion, qr must not be large; values in the range 0.1 
to 0.5 are typical. This constrains the applied r.f. frequency O and hence we obtain 



At the octupole condition, equations H13|) and H17|) give the centre of mass vibrational 
frequencies for motion in the axial and radial directions respectively. Since we want the ions 
to align themselves along the z axis, we require cOr > (and cOr ^ ^^i is desirable). This 
sets a lower limit on the required value of Qac for given (3: 

\UJl J Qr \2TTeo J 

It is useful to examine the way these parameters scale with p, the distance from the origin 
to the nearest electrode surface, and -Bmax, the largest electric field at any electrode surface. 
For a given electrode geometry p sets the distance scale of the electrodes. To keep i?max as 
small as possible, the electrode surfaces should have as large a radius of curvature as possible, 
and therefore their radii should be allowed to increase with p. In other words the complete 
structure of the electrodes is assumed to scale with p. Nevertheless for a given electrode set at 
a given value of -Emax, /? and Qac still have a range of possible values, because it is possible to 
choose a range of values for the ratio between the sizes of the oscillating and the d.c. voltages. 
Typically the variation is such that there is a competition between high (3 and high Qac- 

In order to obtain the main features by a relatively simple analysis, we define two geometric 
factors 7 and /i such that 

forgac^O, P^^^, (19) 

for/3 = 0, Q,, = ^^. (20) 
P 

Thus 7 and p tell us about the maximum values of /3 and Qac available for a given electrode 
geometry, but we keep in mind that it is not possible to have both these maximum values at 
once. 

We can now quantify the desirability of the octupole potential, as compared with other pos- 
sibilities. If the electrode structure is not designed to produce small or zero d.c. quadrupole, 
then it will produce a d.c. quadrupole potential in the xy plane of strength af'^ ~ p'^(3. To 
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achieve radial confinement, we require as given by p7(l to exceed 2af'^q/'m. Hence the 
condition for stable radial confinement is Qac > 4p^/3/(7,- and therefore 

fi > -7. (21) 

Qr 

It is interesting that this condition is scale- independent. Since it is satisfied for large fi or small 
7, it shows that suppressing the d.c. quadrupole in all directions is not necessary for electrode 
structures which naturally produce a strong radial confinement or a weak octupole. However, 
we are interested in making both parameters large, and it will be shown that the structures 
offering high values of both parameters do not satisfy H21|l . therefore the suppression of the 
d.c. quadrupole is useful. 

Suppose the ions each have one electronic unit of charge and let A be the mass number, 
then for a trap at the octupole condition (a = 0), eq. (|13|) gives 

^ir..xTi-., 840 (7gmax)'/^° 

-[MHzlc.— ^^^^^ ■ (22) 

where i?max is in V//im and p is in /im. Note that this mode frequency scales almost linearly 
with and is relatively insensitive to 7 and -Bmax- The radial frequency (eq. HITI) ). on the 
other hand, scales as 

l^r oc (23) 

This difference in the scaling of frequency with trap size, for the radial and axial motion, is 
illustrated in figure|31 The figure shows the frequencies uji (eq. 1)22(1 ') and ujr (eq. (|17l) ') plotted 
as a function of p for two example electrode configurations discussed in section |S1 In order 
that a single graph can show the behaviour for all ion species, the vibrational frequencies are 
shown in units of a basic frequency cuq which depends on the charge/mass ratio of the ion. We 
define ujQ to be the radial secular frequency (eq. IjlTIl ') for an ion in a trap having Qac = 10^ 
Vm-2 and Qr = 0.3. wo/27r takes the value 923 kHz for ^^Ca+, 574 kHz for "iCd+ and 2017 
kHz for 9Be+. 

Each point where the lines cross on figure |31 identifies a distance scale pc- It is found by 
substituting (fT^ and into ((TF|) and solving for p: 

/ 1 \ 1/4 / 1 \ 5/4 / X 5/2 

Pc^'-(-) (^) M (24) 

M \l^J■J 

where 

Lo-36(3g/7reo£;,„ax)'/'. (25) 

The distance scale pc answers the question, "which motional frequency, the radial or the axial, 
will be the more challenging to achieve when we design the electrode structure?" It thus helps 
to establish priorities in the task. Below pc the challenge is to increase the radial frequency, 
while above pc the challenge is to increase the axial frequency. This is because in order to 
maintain a linear ion configuration, below pc the trap must be weakened in the axial direction. 
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P (m) p (m) 

Fig. 3. Centre-of-mass vibrational frequencies for axial motion of two ions in a d.c. octupole 
potential (solid line, u>i) and radial motion of two ions in a linear Paul trap (upper dashed line, 
ujr). The frequency unit luq is defined in the text, (a): electrode arrangement 2.1, (b): electrode 
arrangement 5.1. The frequencies are obtained using values of 7 and /i presented in section 151 
and using i?niax = 10* V/m. The lower two dashed lines are LUrf^ and u!r/5; these permit the 
conditions u!i = uir/2 and wi = LJr/5 to be identified easily. 



by operating it with /? below the maximum available value, and therefore in this regime 
gives the limit on the performance. 

We now show that the value of pc is of order microns in practice. 

In the expression (|24|l the factor {l/6qr)^^^ ~ 0.5 to 2 for a choice of Mathieu Qr parameter 
in the range 0.1 to 0.3. The combination of geometric factors in H24|) has a value in the range 
0.1 to 2 for the structures described below (with the exception of those shown in figures 4 
and 5 which we reject as unsatisfactory owing to their poor radial confinement). Therefore 
we find pc ~(0.3 to 22)Lo for ujr — In practice one would choose LOrji-ox in the range 2 to 
5 in order to keep the ions near the axis without restricting ld\ too much. Increasing Ur/uji 
from 2 to 5 creates an order of magnitude increase in pc- 

In practice we are concerned with singly charged ions, and therefore Lq is a distance scale 
set purely by the maximum electric field allowed in the structure. For example when -Emax = 
10^ V/m (the order of magnitude when field emission may occur Lq ~ 0.15 /zm. This 

field depends mainly on the materials, and partly on the geometry (for example a larger field 
may be allowed across a vacuum gap than across an insulating material surface). However, 
owing to the square root in (|25|l the dependence is not strong. 

Taking qr — 0.3 and combining these observations, it is seen that for a variety of structures, 
Pc < 10 pm for uJr = 5uji at high electric field, and pc < 2 pm for ujr — 2cji and more modest 
fields. This implies that for ion traps designed for quantum computing experiments, the 
more demanding problem is currently to obtain a large d.c. octupole moment, but the radial 
confinement may become the more demanding problem in the future as distance scales are 
reduced. The values of p where the lines cross in figure |31 illustrate these points. 

Another important distance related to Lq is the distance between the ions when the 
electrodes are close to the octupole configuration. Using 10 and p9|l we obtain 
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Therefore d p when p 3> Lq/j^^'^, which is usuaUy the case. For example, putting 7 — 0.02 
(c.f. table 2 below) and Lq ~ 0.2 fim, d <^ p when p ^ 16 nm. For electrodes of micron 
scale or above, therefore, the available confinement is such that when the electric potential is 
just at the point of producing a potential hill between the ions, the distance between them is 
much smaller than the distance to the electrodes. This implies that the Taylor expansion we 
have adopted is valid. 

3. 1 Scaling of heating effects 

The main reason why one would not make a quantum computer faster simply by fabricating 
the traps at a smaller distance scale is the fact that heating problems will be exacerbated. 
Heating rates observed in various studies scaled as p~'^ ^1 and the scaling expected for 
heating caused by Johnson noise is p^^ . We will treat some general scaling law p^^ , and 
we assume this heating depends on the electrode structure only through the distance scale 
p. The time to separate a pair of ions scales as l/wi, therefore the number of phonons' 
worth of heating per split time. An, is proportional to 1/(^1/9'^). This gives a good guide to 
the amount of cooling which will be needed to obtain a steady average temperature in the 
computer. Using equation (|22|l we obtain 

-fc+9/10 

^^^-rprsT- (27) 

For a given electrode geometry, 7 is fixed, hence this scaling law governs the allowed distance 
scale p. If we take An as a fixed quantity, then p oc ^3/(9- 10'^). Substituting this behaviour 
into equation H22|) we find the following scaling for the trap frequency in terms of 7: 

c^i cx 73'=/(io'=-9) (28) 

For example, for fc = 2 the scaling is oji cx 7"-545^ for fc = 4 the scaling is wi cx 7"-387^ When 
the radial rather than the axial frequency is the harder to achieve, then a similar analysis 
may be used to relate Ur to p. 

3.2 Study of a non-octupole trap structure in current use 

An example of an ion trap in current use which has been used to separate ions is that operated 
at NIST Boulder for teleportation and other experiments K is a two- layer trap with 5 
d.c. electrodes in two opposing quadrants, and the r.f. electrodes in the other quadrants. 
With z along the trap axis, the ends of the electrodes lie at x = a, y = a where a — 100 /im. 
It is not possible to produce an octupole using such an electrode configuration (this requires 
7 d.c. electrodes, see section 1^^ . 

We studied this structure by modelling it and solving Laplace's equation numerically. The 
calculations were performed using Simion and CPO [53]. Simion uses a finite difference 
method to find the potential of a regular 3D array of points. The array is iteratively relaxed 
so as to satisfy Laplace's equation and the boundary conditions, until the fractional changes 
per iteration are all below 5 x 10"^. However, from simulations of analytically solvable cases, 
we found that the error in electric fields close to spherical and cylindrical electrodes could be 
as high as 10%. CPO uses a surface charge method, and we found it to be considerably faster 
and also much more precise. All the numerical results presented in this paper were obtained 
using CPO. 
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To calculate a value for the parameter 7 = 13 / E-a^s^, we numerically solve Laplace's 
equation for the given electrode set-up, and then extract /3 — {\/2A)d'^V/dz^ evaluated at the 
origin, and .Emax by examining the gradient of V near the surfaces of all the electrodes. 

In order to find /z, calculations were performed with the r.f. electrodes at a voltage of 
10^ V which is large compared to the d.c. voltages of ~ 1 V used for the d.c. electrodes. 
The value of the radial confinement = Qac was extracted, along with the value of the 
maximum field E!^-^^^ at the surface of the electrodes. The radial geometrical parameter /i was 
then calculated using equation H2Q|I . 

For the model of the NIST trap, first the r.f. electrodes were grounded and the voltages 
on the d.c. electrodes were adjusted to satisfy d'^V/dxdy = d'^V/dydx = d'^V/dz'^ — 0. We 
then found that d^V/dy"^ = -d'^V/dx^ = 0.02 Va'^, (i = 0.01 Va"* and £;,„ax = H-l Va"!. 
The NIST trap has p = 1.38a, hence we deduce 7 = 0.0024. The same procedure was carried 
out with d^V/dx^ = d^V/dy^ — d^V/dz^ — but the remaining (undesirable) quadrupole 
terms were found to be 10 times larger. In each case we then raised the r.f. electrodes to a 
high voltage in order to find /j,, obtaining the value 0.11 in both cases. 

We observe that the d.c. quadrupole which exists in the x ~ y plane when d'^V/dz'^ ~ 
in this trap has a strength very close to i.e. it is just as large as we would expect it 
to be from general scaling considerations when the structure is not designed to eliminate it. 
The values of 7 and /i satisfy condition (|21(l for > 0.006, therefore in this trap the radial 
confinement from the a.c. potential can overcome the d.c. anticonfinement before the r.f. 
voltage is restricted by electrical breakdown. This can be regarded as owing to a low value 
for 7. However, the presence of the d.c. quadrupole remains unwelcome, because it results in 
a large increase in the required r.f. voltages. 

To estimate the maximum axial and radial frequencies which could be achieved at the 
point of ion separation, we use qr — 0.3 and assume that the maximum field allowed at the 
surface of the electrodes is 10* Vm~^. Then the voltages on the d.c. electrodes are 67, - 
153, 509 V on the central, intermediate and outer electrodes respectively. For tOr > 3^^ the 
r.f. voltage amplitude \s Vrf = 557 V. In the case of beryllium ions, these conditions lead 
to motional frequencies (w^: , "^y , ) = 27r x (7,30,2) MHz. The large difference between lo^ 
and ujy shows that most of the voltage on the r.f. electrodes is being used to overcome the 
d.c. anticonfinement along x. If the d.c. quadrupole were eliminated, then one would obtain 
uj^I'l-K ^LOyjl-K both equal to 20 MHz. In that case one would have the option of reducing the 
r.f. voltage by an order of magnitude. 

4 Realising octupoles 

In this section we discuss general methods to obtain an octupole moment at the origin of 
coordinates. This acts as a guide in the design of the electrodes. 

Three examples of octupole potentials are the axially symmetric octupole: 




(29) 



an octupole with cubic symmetry: 



T/ 4 I 4 I 4 o/' 2 2 I 2 2 I 2 2\ 

K = x + y + z — S[x y + X z + y z ) 



(30) 
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and the 2-diinensional octupole in the x-z plane: 

V = z^ + x'^ -6z^x^. (31) 

The axiaUy symmetric octupole can be produced by a set of 2 (+ve) end-cap and 3 ring (2 
-ve, central one +ve) electrodes, all shaped to follow contours of V. The 2D octupole can 
be produced by a set of 8 (4 +ve, 4 -ve) rod electrodes parallel to the y axis and located at 
the corners of an octagon. These two examples of pure octupoles give a useful pointer to the 
required distribution of image charge. 

The methods and image charge distributions listed below do not produce pure octupoles, 
but have an octupole as the leading term in the multipole expansion of their potential near 
the origin. 

4-1 Use of symmetry 

To produce an octupole at a point r, 19 constraints have to be satisfied. These are that the 
first, second and third derivatives of the potential at this point in three orthogonal directions 
should vanish. Many of these can be satisfied by the introduction of symmetry to the charge 
distribution producing the potential. This can take the form of cither rotation symmetry 
about an axis passing through r, or reflection symmetry in a plane containing the point r. 

A single reflection in a plane containing the point r produces zero odd-order moments in 
a direction perpendicular to the plane, and also forces the mixed second derivatives such as 
d^V/dxidxj to vanish, where Xi lies within and xj lies perpendicular to the plane of reflection. 
This leaves 12 constraints to be satisfled. 

Two reflections in orthogonal planes containing r produces zero odd-order moments in the 
directions normal to the two planes, and all mixed second derivatives vanish. This leaves 7 
constraints to be satisfied. 

Two-fold rotation symmetry about an axis passing through r causes the odd-order deriva- 
tives of the potential in all directions perpendicular to the axis to vanish, reducing the number 
of constraints to 9. 

Four-fold rotation symmetry about an axis passing through r makes the second derivatives 
in directions perpendicular to the axis equal, hence the number of constraints is reduced to 
5. 

Combinations of the above can be used to reduce the number of constraints further. 
Reflection symmetry in three orthogonal planes all containing r reduces the constraints to 3, 
on the non-mixed scc;ond derivatives d'^V/dxf. Combining a two- fold rotation about an axis 
passing through r with a reflection in a plane normal to the axis and containing r reduces 
the constraints to 4 (the non- mixed second derivatives and the derivative d^V / dxidxj, where 
Xi and Xj are orthogonal directions perpendicular to the rotation axis). A four-fold rotation 
symmetry about an axis passing through r combined with a reflection in a plane normal to 
the axis and containing r leaves 2 constraints (on the three non-mixed second derivatives, two 
of which are equal) . 

Finally, Laplace's equation reduces the degrees of freedom of the non-mixed second deriva- 
tives d'^V/dx'l by one. 

4-2 Image charge constructions 

(A). An octupole with cubic symmetry is produced by any arrangement of charge having 



14 Electrode Configurations 



cubic symmetry. Proof: the symmetries under reflections in three orthogonal planes imply 
the odd-order multipole moments vanish, and so do mixed second derivatives of V such as 
d^V/dxdy; the symmetries under rotation through 90° imply the coefficients of x^, and 
in a Taylor expansion of the potential at the centre of symmetry are all the same, and 
therefore by Laplace's equation they all vanish. 

(B) . An octupole with cylindrical symmetry is produced by any arrangement of charge 
having cylindrical symmetry, reflection symmetry in z = 0, and placed such that d'^V/dz'^ = 
(where z is the symmetry axis). Proof: the second derivatives w.r.t. x and y must be equal 
by symmetry and therefore by Laplace's equation they also vanish. 

(C) . An octupole is produced by any arrangement of charge having four-fold rotational 
symmetry in the x — z plane, reflection symmetry in y = 0, and arranged so that d^V/dy^ = 0. 
The reasoning is the same as in the previous case. The resulting octupole potential is a 
combination of the axially symmetric case H29|l and the 2-dimensional case oriented such 
that y is the common axis. 

(D) . An octupole can be produced, starting with any arrangement of charge p{x,y,z), 
by the following recipe. First ensure reflection symmetries in order to cancel odd multipole 
moments, by forming 

111 

^ ^ E E E pii-^y^^ i-^yy^ (-i)'^)- (32) 

i=0 j=0 k=0 

Next, form 

p' = p{x,y,z) - p{x/f,y/f,z/ f). (33) 

where / ^ 1 is a numerical factor, p' produces an octupole at the centre of symmetry. This 
can be proved by evaluation of the second derivatives of V{r) oc J p'{r')/\r — r'jdr at the 
origin. A proof giving further physical insight is as follows, pix/ f,y / f, z/ f) produces the 
same functional form of V{x, y, z) as p{x, y, z) but on a distance scale expanded by a factor 
/. Modelling p as a set of point charges, each charge in the rescaled distribution is now 
further from the origin by a factor /, and is increased in size by a factor f^. Since the second 
derivatives of V due to a point charge go as (distance)"'^, it follows that the second derivatives 
of the new potential are equal and opposite to those of the original one, QED. The fourth 
derivatives scale as (distance)"^, so 

where V' (V) is the potential due to p' [p) respectively. 

A natural extension of this method can produce higher order multipole configurations. 

(E) . An octupole can be formed by taking almost any charge distribution p(x,y,z), and 
displacing and reflecting it in three dimensions. A general displacement to p{x — xo, y — yo,z~ 
zq) has three free parameters (xcj/Oi^o)- One can be absorbed into an overall scale factor, 
leaving two. This is just enough free parameters to allow the two constraints 



dx"^ 



(0,0,0) = 
(0,0,0) = (35) 
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to be satisfied, where V is the potential due to p(x — xo,y — i/q, z — zq), except in rare cases 
when the form of p leads to no solution. The nulling of d^V/dz^ follows from Laplace's 
equation. Finally, introduce reflection symmetry as in (jSU, which causes odd-order moments 
and mixed second derivatives {d^V/dxdy etc.) to vanish. 

Another way to understand this construction is to argue that most charge distributions 
produce a potential with zero second derivatives with respect to fixed orthogonal axes at some 
point in space; it suffices to place the origin at such a point and then introduce reflection 
symmetry. For example, if p corresponds to a single point charge, this construction leads to 
point charges on the corners of a cube. More generally if p{x, y, z) lies in the x-z plane, this 
construction leads to a distribution of charge in two parallel planes at y = ±2/o- 

(F). All the above constructions involve constraints on the locations of the image charges 
in addition to the basic assumption of three-fold reflection symmetry. Our next construction 
does not. An octupole can always be formed (with rare exceptions arising in pathalogical 
cases) by a charge distribution p with reflection symmetry as in ^A'l\ . based on 



where pi{x,y,z) (i — 0,1,2) are three different charge distributions located anywhere in 
the positive octant, and /i and /2 are parameters whose values are given by solving the 
simultaneous equations 



(evaluated at the origin), in which Vi are the potentials due to p,;. This construction uses 
the same property invoked in the previous one, namely that after cancellation of odd-order 
moments there are only two constraints Ij^Sfl we need to satisfy. This implies that we only 
need two free parameters, and these are provided by /i and /2. 

(G). When we abandon a reflection symmetry, the octupole can be regained by inserting 
further charges. For example, suppose we only have reflection symmetry in the x~y and y-z 
planes. This is the case, for example, when all the charges are located in or on a single 
substrate, with the octupole centred above the top surface of the substrate. Then in addition 
to we have 4 further constraints, on dV/dy, d^V/dy^, d^V/dx^dy, d^V/dz'^dy. One 
or more of these can be satisfied by a careful placement of a single charge distribution. 
Alternatively, charges are placed at any convenient position, and then their magnitudes are 
adjusted, similar to construction (F). To obtain the desired behaviour of the trapping centres, 
we don't need to insist on an octupole, however. Of these further constraints, only dV/dy = 
is strictly necessary. 

4-3 Application to electrode design 

The distributions of charge discussed above are to be realised on the surfaces of a set of 
conducting electrodes at fixed voltages (or else they are image charge distributions which lead 
to the same equipotential surfaces). We wish to avoid complicated electrode shapes so we 
restrict ourselves to electrodes approximating to simple combinations of lines, rings, or sheets. 

In order to satisfy N constraints on derivatives of V{x, y, z) merely by adjusting voltages 
on electrodes, it is necessary to have -t- 2 electrodes, since one defines a voltage zero, and 



p{x, y, z) = pn + fipi + /2/O2 



(36) 




(37) 
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the next one sets an overall scale factor. In order to produce radial confinement, a set of 
r.f. electrodes is introduced around (or in between) the d.c. ones which create the octupole. 
These r.f. electrodes also serve to define the d.c. zero. Thus, in construction (F) we require 
one r.f. voltage at d.c. ground, and 3 adjustable d.c. voltages; in construction (G) we require 
one r.f. voltage at d.c. ground, and 7 adjustable d.c. voltages. 

Most of the constructions presented in the previous section can be ruled out, as follows. 

The case of cubic symmetry (A) involves 8 electrodes, but to preserve precise symmetry 
these would have to lie along the edges or the extended major diagonals of a cube, which is 
difficult to fabricate, and a linear r.f. electrode set would also break the symmetry. 

Case (B) can be realised with just two ring-shaped electrodes. This involves the fewest 
electrodes. For a Penning trap one would choose the magnetic field direction as the axis of 
symmetry of the rings. For a Paul trap the need to provide an additional r.f. quadrupole 
producing radial confinement in the x — y plane suggests that there will be further electrodes 
parallel to the z axis. This implies the fabrication is more straightforward if the rings are 
placed around another axis, which we take to be the y axis, so that all the electrodes lie in 
a set of parallel planes. The r.f. electrodes then destroy the circular symmetry. To regain a 
simple design, the rings can be replaced with rectangular loops, c.f. construction (C). 

The other cases can all be fabricated in layers. The method of (case (D)) involves 
up to 16 electrodes in total or 8 if they all lie in a plane; the method by displacement (case 

(E) ) involves 8 electrodes in two parallel planes, the method of H36|l (case (F)) involves up to 
24 electrodes in general, or up to 12 if they all lie in a plane, which we take to be the x — z 
plane. This can be reduced to 10 if one pair lies on the x axis. We avoid placing electrodes 
along the z axis, because we want to be able to bring ions into and out of the structure by 
moving them along this axis. 

The minimal number of electrodes for a perfectly constructed system is not the only 
consideration, however. Manufacturing imprecision, patch potentials, and the desire to adjust 
the trap separation all imply that further electrodes will be needed to trim the system. The 
electrode structures which rely on precise construction to produce an octupole (cases (A)-(E)) 
require more of these extra electrodes than those which depend on variable voltages (cases 

(F) and (G)). 

Screening effects arise when we consider conducting electrodes rather than merely charge 
distributions. The innermost electrodes tend to act like a Faraday cage to modify and screen 
the effects of the outer ones. Let P be the centre of the region in which we require a certain 
potential form (here, P is at the origin). A rough rule of thumb is that if there is not an 
unobscured line of sight from P to a given electrode, then the influence of that electrode 
on the potential at or near P is greatly modified by the obstructing electrode. However, 
the outer charge distribution in (|33l) is, by construction, exactly on an extended line from 
the origin to the inner charge distribution. Therefore an electrode set built with the aim of 
realising the recipe (D) suffers from screening effects and it is difficult to obtain a surface 
charge distribution of the form Ip^ . 

In the remainder of this paper we will concentrate on the case (F) of 10 electrodes all in 
the x-z plane, of which one pair lies along the x axis, the case (E) of 8 or more electrodes 
in two planes at y = ±yo, and the case (G) of an octupole situated above all the electrodes. 
For the sake of completeness, the appendix lists some other examples of octupoles based on 
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simple distributions of point, line and ring charges. 

5 Optimization of some example structures 

In this section we study various electrode arrangements and obtain the geometric factors 
7 (eq. (23) and ^ (eq. 1(201)) in each case. We optimize the most promising designs by 
adjusting electrode positions and relative sizes in order to maximise both factors. 

In all the configurations to be discussed, we choose the z axis along the line where the 
twin traps are situated, the x — z plane is parallel to the planes containing sets of electrodes, 
and hence these planes are separated along the y-direction. 

The geometric factors 7 and were obtained by the numerical method described in section 
13.21 To obtain ^ the r.f. electrodes were assigned a large voltage and the second derivative 
of V{x,y,z) was calculated along the principal axes of the quadrupole perpendicular to the 
z axis. In most cases these were the x and y directions (giving and ^y). Where this is 
not the case the radial geometrical parameters are denoted by and ^.yi. Where and pLy 
differ, a component of the oscillating quadrupole potential lies along z, hence trapped ions 
would experience axial micromotion and the associated pseudopotential would contribute to 
their confinement along the z direction, as discussed at the end of section 12.21 

Since the numerical solution of Laplace's equation can be slow, we first gained some general 
insights by analytic methods. We roughly modelled the electrodes as a set of line charges, 
and wrote down the complete potential function. This permits 7 to be found as a function 
of various parameters, such as the line charges and their positions, and it can be maximised 
analytically. Results from such analyses gave a starting point for the electrode structures 
used in the numerical studies. 

In order to know what value of 7 and /i one might hope for, we numerically calculated 
these parameters for two simple cases. These can be regarded as a 'standard' against which 
the results for our various electrode structures can be measured. As a standard for ^, we used 
four cylindrical electrodes of radius 0.41254a with axes aligned along the z axis, and centred 
at the corners of a square of side 2a. This arrangement gave ji = 0.20. Electrodes of radius 
0.3a centred at the same positions gave /i = 0.15. 

To obtain a standard value for 7, an electrode structure loosely based on the equipotentials 
of the axially symmetric octupole in equation 1291 was used. This consisted of three toroidal 
electrodes centred on the z axis at z = and z = ±0.643a, with two spherical "end-cap" 
electrodes of radius 0.35a at z = ±a. The cross-section radii of the torii tubes and the radii of 
the centres of the tubes were 0.2a and a respectively for the z = torus and 0.2a and 0.766a 
respectively for the z = ±0.643a torii. This configuration gave 7 = 0.133. 

5.1 8 d.c. electrodes in 2 planes, 2 r.f. electrodes 

Let us first consider structures consisting of rod d.c. electrodes placed at a distance p 
above and below the plane containing the r.f. electrodes. This "sandwich" arrangement is 
shown in figure 0] a). It is based on construction (E) of section IT!^ All the d.c. electrodes 
are at the same voltage. 

Plan views of the two d.c. electrode structures tried are shown in figures^ and[Sl We let 
a set the overall distance scale, and then there are 4 parameters (6, d, p, r) and 2 constraints. 
The octupole configuration was found by fixing r and b, and then adjusting p and d. 
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Fig. 4. Views of the electrode structure with 4 d.c. electrodes in each of two planes, and 2 r.f. 
electrodes, a) cross-section in the x — y plane showing the sandwich structure. The r.f. electrodes 
are filled grey, b) View looking down on the x — z plane from above. 



The equipotential along the r.f. electrodes tends to prevent an octupole being produced 
at the origin. We found that it was necessary to restrict the influence of this equipotential by 
keeping p small. Since p must be larger than the electrode diameter this means that r must 
not be too large. 




Fig. 5. A view of the electrode structure with 45 degree rod electrodes having hemispherical ends 
centred at (ita, itp,±d). The r.f. electrodes have thickness 2r, and lie in the y = plane. The 
end of these electrodes are hemicylinders of radius r and centered on a; = ±6, y = 0. The grey fill 
represents the r.f. electrodes. 

Results of the calculations for both structures are presented in table ^ The two cases 
produce similar values of 7 and /i, those (labelled 1.1) for the structure shown in fig. being 
somewhat higher than those (labelled 1.3) for the structure shown in fig. [3 It was found that 
the r.f. electrodes had to be placed well back (a large value for b) in order to get a high value 
for 7. 

Results 1.1 and 1.2 show that the value of /i decreases as b is reduced. This is because p 
and d also decrease, with the result that the influence of the r.f. electrodes is screened by the 
d.c. electrodes. 
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The geometry shown in figure ^provides greater screening of the r.f. electrode by the d.c. 
electrodes, compared to that shown in figure El and this is probably the main reason why it 
provides the higher 7. 

Table 1. Values of the distance parameters r,b,p,d which give an octupole for three example 
electrode structures, with the distance to an electrode surface p and the calculated geometric 
factors 7 and fi in each case. Results 1.1 and 1.2 are for electrode configurations of the form shown 
in figure |4] Result 1.3 is for an electrode configuration like that shown in figure l5l 
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5.2 10 d.c. electrodes in a plane, 4 f-f- electrodes 

This structure consists of 10 d.c. electrodes in the y = Q plane, and 4 r.f. electrodes in two 
planes parallel to this. It has reflection symmetry 'm x — y, x — z and and y — z planes. The 
d.c. electrodes are cylindrical, with hemispherical ends centred at a; = ±a. Figure shows a 
cross section in the x — z plane for positive x. 
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Fig. 6. 10 d.c. electrodes, 2 r.f. electrodes: A cross section view of the d.c. electrodes in the x — z 
plane with y = for positive x. The trap centre is at the intersection of the x and z axes. The 
electrode structure is repeated for negative x by reflection in the z — y plane. 

Each r.f. electrode is effectively a semi-infinite plane, with its edge parallel to the z axis. 
The y position of the r.f. electrodes is zkp, and the x-position is zL{a + t), hence t gives 
the retraction of the r.f. electrodes. Where p = 1.03 the r.f. electrode had thickness a/4, 
otherwise the thickness was negligible. For the finite thickness cases, the edge of the electrode 
is hemicylindrical, and p and t give the position of the axis of the hemicylinder. An electrode 
structure with p = 1.03a and t = a is illustrated in figure 6. 

The voltage on electrode 3 (V3) was set to 1.4V in all the calculations in this section. 
With a fixed geometry our free parameters are the voltages on electrodes 1 and 2 (Vi and V2 
respectively). These were adjusted to fulfill the octupole condition, using a two-dimensional 
Newton-Raphson method in order to find the solution with few trials. Then the values of 
Emax and 7 were calculated. Adjustments were then made to /, d and c?3, and the calculation 
repeated. A selection of positions and their corresponding values of 7 are shown in table [21 
The nearest electrode surface distance p = 0.825a in all cases. 
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Fig. 7. A perspective view of the electrodes for r.f. electrodes positioned at p = 1.03a and t = a. 
The 10 cylindrical electrodes have d.c. voltages applied to them. r.f. voltages are applied to the 
planar electrodes. Also shown are insulating spacers which hold the electrodes in place. 

The initial choice of electrode diameters was di = 0.35a, d2 — 0.825a, = 1.3a, based 
on the indications from our initial line charge calculations. Using these diameters, CPO 
calculations were carried out with no r.f. electrodes for several values of / and d. The 
optimal values were found to be / = 3.25a and d = 0.825a. 

We then introduced the r.f. electrodes. With the r.f. electrodes far away from the d.c. 
electrodes {p = 2a and 4a) the maximum field is found between the d.c. electrodes, hence 
varying / makes little difference to the value of 7. This is because an increase (decrease) in 
j3 is countered by an increase (decrease) in i?max- Cases 2.1 and 2.2 in table |21 illustrate this. 
For large d, the maximum field is found between electrodes 2 and 3, hence as d is reduced, 
Emax decreases. For small d, the maximum field is found between electrodes 1 and 2, and 
then as d is reduced, E^^^y^ increases. This behaviour is shown in figure |H1 

Table 2. Results from CPO calculations. Results 2.1-2.6 are for the structures described in section 
|01(c-f- figures|nil3, results 3.1-3.3 are for the inverted structures described in section 15.31 All 
the cases listed have d = 0.825a, di = 0.35a, d-i = 0.825a, V^f = and V3 = 1.4V. The table 
illustrates effects of changing the values of /, d^, p and t. 
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When the r.f. electrodes are close to the d.c. electrodes {p — 1.03a), the maximum field 
is found between the d.c. electrode 3 and the r.f. electrodes. This field can be reduced by 
reducing ds. (3 increases as / is reduced, hence 7 can also be increased by reducing / until 
the maximum field is found between the d.c. electrodes. Figure El shows this behaviour for 
structures with p = 1.03a and i = 0, and the net effect is illustrated by results 2.5, 2.6 in 
table El 

We also examined the effect of making electrode 3 extend further along the z direction. It 
was replaced by a planar electrode of thickness d^ with edges rounded off by hemicylinders, 
located such that the edge nearest the origin was situated at the same place as the cylindrical 
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Fig. 8. The effect on /3, Emax and 7 of adjusting d for result 2.1 in table |5] 100/3 is represented 
by +, and Emax by unfilled circles. If the largest field lies between the first and second electrodes, 
then reducing d increases both /3 and -Emax, and 7 drops. Where this is not the case, the increase 
in /3 causes an increase in 7. 




b) 




Fig. 9. The effect on /3, Emax and 7 of adjusting /, with d = 0.825a, di = 0.35a, ^2 = 0.825a, 
da = 0.825a, r.f. electrodes at p = 1.03a, 4 = 0. 100/3 is represented by +, and Emax by circles. 
Reducing / has the effect of increasing /3. At larger /, the maximum field is found between the r.f. 
electrodes and the d.c. electrodes, hence we can increase 7 by decreasing /. After reaching the 
point where the largest field is found between the d.c. electrodes, 7 is much reduced for further 
reduction of /. 
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electrode in the other results. This had little effect on 7 or ^. 

Finally, we examined an example in which the d.c. electrodes are extended along the z 
direction into planar shapes. This will reduce 7, but it also reduces the z-dependence of the 
a.c. part of the potential, therefore giving less micromotion along z. Each electrode had 
two flat surfaces parallel to the x — z plane, and sides in the shape of hemicylinders having 
the same radius of curvature as before. The electrode centres were placed as in result 2.6 
and they were extended along the z direction until the gaps between them were 0.1a. This 
arrangement gave 7 = 2 x 10"'^, i.e. a factor 10 reduction compared to result 2.6. 

5.3 20 d.c. electrodes in two planes, 2 r.f. electrodes 

This arrangement is like the previous one, but with the d.c. and r.f. electrode positions 
swapped. It consists of 2 planar r.f. electrodes of thickness 0.5a lying in the y = plane. 
These are bounded by hemicylindrical edges centred at a; = ±(a + t). The d.c. electrodes 
lie in two planes at y = ±p, each containing two sets of five d.c. electrodes (with /, d and 
electrode diameters the same as result 2.1 in table EJ- Results for these arrangements are 
shown in table [51 as results 3.1-3.3. The values of 7 and fi are similar to those in section 
Some improvement in 7 would be expected through optimisation of the relative positions of 
the d.c. electrodes. This has not been carried out in obtaining results 3.1-3.3. 

5.4 Low Aspect Ratio Electrode Configurations 

Let w be the total thickness, i.e. depth in y, of the electrode structure (c.f. figure [Tn|. We 
define the aspect ratio g = w/p. The three-layer trap designs considered in sections l5. II 15.21 
and 15. 31 have g > 1. Though microfabrication techniques are available which could be used to 
construct traps with p < 100 pm and g > 1 |26j . many methods are limited to a maximum 
structure thickness of around 20/im |16l 12 7j . This means that traps may have a low aspect 
ratio, g < I- 

We studied structures similar to those of section but now with .g < 1. To speed the 
calculation process we used cuboid shapes for the electrodes (see figures [TUl and ITTll . The d.c. 
electrodes had width in z of di = ^2 = da = w, and were separated by gaps of w/2 (hence 
d — 3w/2, f = Sw). The thickness of the electrodes in y was w/5. With a cuboid shape the 
electric field diverges near the electrode edges. In order to calculate 7 and p we need to make 
some reasonable assumption to handle this. Our aim was to use the calculation to indicate 
approximately what would happen if in fact the electrode edges were smoothed. To this end 
we calculated electric field values at (x, y, z) = (p, 0, 0), {p+w/^, 0, di/2), (/9+t + w/5, iy/10, 0) 
and at similar places around the other electrodes, and took £',nax to be the largest of these. 

Results for structures with g = 0.4 and g = 0.2 are shown in figure [T^ We found that 
retracting the r.f. electrodes (i.e. using t > Q) gave a useful increase in both 7 and p. As the 
aspect ratio decreases, 7 falls rapidly (and p less rapidly). 

5.5 Two-Layer Electrode Structures. 

The most natural configuration for producing a radial r.f. quadrupole in the x — y plane 
consists of electrodes placed at four corners of a square. Equal r.f. voltages are applied 
to diagonally opposite electrodes. To produce an octupole by such an arrangement, one or 
more of these electrodes must be segmented. Choosing the structure to have 2-fold rotational 
symmetry about the z-axis, and reflection symmetry in the plane z = 0, there remain 3 
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Fig. 10. A cross-section view of a low aspect ratio electrode structure for positive x showing 
the definition of the parameters ui and t. The r.f. electrodes are represented by grey fill. The 
electrodes for negative x are obtained by a reflection in the x = plane. 



Fig. 11. Perspective view of a low aspect ratio electrode structure with g = 0.4 and t/ p = 0.6. 
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Fig. 12. Values of a) 7, b) for low aspect-ratio traps, as a function of the distance the r.f. 
electrodes are retracted with respect to the d.c. electrodes. Values for g = 0.4 structures are 
marked by *, and values for g = 0.2 are marked with a -f-. For a given g, the lower values of fi are 
p,y, the upper values are p,x- 
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constraints. We chose to satisfy the constraints by adjusting vohages rather than by accurately 
placing electrodes, and therefore two diagonally opposite electrodes must be split into seven 
segments. 



a) 
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Fig. 13. The segmentation of the d.c. electrodes used in the two-layer electrode structures for a) 
cuboid electrodes and b) cylindrical electrodes. For cases where cuboidal electrodes were used, 
the electrode widths di to d4 were equal, and the separation between electrodes was di/2. For 
structures with cylindrical electrodes, d2 = ds = d4 = 2.36di. 



The d.c. electrode segmentations used in the results presented below are shown in figure 
[T^ The cuboid electrode structures consist of electrodes of equal width {di =^2 = ^3 = ^4), 
equally spaced by c?i/2 (hence d — 3/2di, f — 3d and h = 9/ 2d). The cylindrical electrode 
structures were made up of cylinders parallel to the x axis with hemispherical ends. The 
diameters of the cylinders were d2 = d^ = d^ ~ 2.36(ii. The axes of the cylinders were 
positioned at d = 2.36di, / = 5.43di and h ~ %A?>di. The hemispherical ends of the cylinders 
were centred on x = a. 

The r.f. electrodes consist of extended planes of thickness Wrf. For structures with cuboid 
electrodes, the end of the r.f. electrodes were placed at x = ±(a + t). For structures with 
cylindrical d.c. electrodes, the r.f. electrodes had hemicylindrical ends of radius Wrf/2 with 
axes parallel to the z axis aX x — ±{a -\- 1) . 

The voltage on d.c. electrode 4 was set to 1.4V, and the r.f. electrode voltage was set to 
OV. Numerical calculations were performed for a range of values of c?i, t and g, as defined in 
figures [T^ and 1141 Example results of these calculations are presented in tables |2I and 0] 

We found that the width di of the central electrode had a strong effect on 7 and a smaller 
effect on /i (see figure [T51 and c.f results {4.5,4.6}, {4.9,4.10}). The optimum value of di for 
producing a d.c. octupole was found to be around 0.4p, for all values of g. This is comparable 
with the value used for the three layer electrode structures in section [5.21 which was obtained 
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a) 



b) 



c) 



Fig. 14. Profiles of the two-layer electrode structures for negative x viewed in the x — y plane. The 
r.f. electrodes are represented by grey fill, a) is the structure used to obtain results 4.1 to 4.12 in 
table|21 b) was used for results 4.13 to 4.15, and c) for results 4.16 and 4.17. The aspect-ratio of 
the trap was defined hy g = w / p. The structure is repeated for positive x by a rotation of 180° 
about the z axis. 



by examining the equipotentials of line-charge calculations. The monotonic increase in /i with 
di may be ascribed to an increased influence of all of the d.c. electrodes. 
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Fig. 15. The effect on a) 7 and b) fi of adjusting the width of the central electrode for a two-layer 
electrode structure having g = 2 and i = 0. fj,^/ is represented by unfilled circles and fiyi by +. 



Figure [TCI shows the effect on 7 and fi of changing the aspect ratio of the trap. For g < 1, 
both parameters increase with g. At higher g, ^ continues to increase with g while 7 decreases, 
presumably due to the reduced shielding of the r.f. electrodes by the d.c. electrodes. 

We found that, for g = 2, retraction of the r.f. electrodes had no significant effect on 7, 
but reduced /i somewhat (results {4.1-4.3}). When g = 0.4 retraction allows a small increase 
in 7 with a drop in (results 4.7,4.8). When g = 0.24, increasing t from to p/2 increased 
7 by 40% with little effect on ^ (results 4.10,4.11). 

We examined the structure shown in figure [Hh. with the d.c. electrodes placed in the 
y — plane, with a view to finding out whether it would produce a higher 7. However, 
comparison of results 4.17 and 4.15 suggest that it has the reverse effect. Since fi decreases 
also there is nothing to be gained by adopting this structure. 
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Fig. 16. The effect on a) 7 and b) fi of adjusting the aspect ratio g = w/p for a two layer 
electrode structure. In b), the values of fj, plotted are those along the x' (represented by +) and 
y' (represented by o) axes of the quadrupole. The values of di for these cases are similar but not 
identical (0.39, 0.47, 0.4 and 0.4 p from left to right on the plots). 



5.6 Planar or near-planar arrangments 

We considered cases where all the electrodes of the twin-trap system lie in a single plane 
(figure [T7|) . and a simple 'railway track' arrangement in which rod-shaped r.f. electrodes run 
orthogonal to a set of simple rod-shaped d.c. electrodes (figure [TS|) . The simplicity of such 
designs gives obvious advantages from the point of view of microfabrication, but we find that 
the penalty in terms of reduced factors 7 and fi is severe. These results are presented in table 

E 

The voltages on these structures are to be adjusted to produce a trapping region centred 
outside the electrode structure. Therefore they lack one reflection symmetry, and there are 
6 constraints if we require an octupolc (construction (G) of section I4.2|l . To satisfy all these 
constraints by adjusting voltages would require 7 d.c. voltages (in addition to the d.c. ground 
set by the r.f. electrodes) and therefore a large number of electrodes. In order to reduce the 
number of electrodes, we chose to examine structures such that the hexupole moment may 
be small but non-zero. We adjusted the 5 voltages labelled A-E in figures [T7I and IT^ so as to 
cancel the quadrupole terms and dV/dy and d'^V/dy^, but we did not constrain d^V / dx^dy 
or d^V/dz^dy. 

To find the condition where the d.c. hexupole/octupole is centred at the same height as 
the r.f. quadrupole, we proceeded as follows. First the r.f. electrodes were set to zero, and 
the d.c. voltages were adjusted to produce a hexupole/octupole at a few different heights. 
Next, the r.f. electrodes were allowed to go to ±1 and the quadrupole position identified. We 
could then make a good first guess of the right d.c. voltages to get the desired coincidence, 
and the solution was found iteratively. 

6 Manufacturing imprecision and stray electric fields 

Throughout the above we relied on symmetries in the electrode structure in order to 
achieve the desired electric potential function. It is necessary to check whether there might 
be a high sensitivity to manufacturing imprecision and patch potentials which will in practice 
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Fig. 17. The planar electrode arrangement (result 5.1). The cylindrical electrodes all lie in the 
same plane and have the same radius of 0.3a. The ends of the electrodes parallel to the x axis are 
hemispherical. 
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Fig. 18. Railway track electrode arrangement (result 5.2). The electrodes parallel to the z axis 
are situated above those parallel to the x axis by a distance of a. All electrodes are cylindrical, of 
radius 0.3a. 
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break such symmetries. To this end we studied an example case: the electrode structure 
leading to result 2.2 in tabled c.f. figureEl 

In order to illustrate manufacturing imprecision, the electrode which has its spherical end 
centred at (a, 0, d) was displaced by 0.05p in the z direction. The d.c. electrode voltages were 
then adjusted to satisfy d'^V/dx'^ = d'^V/dy'^ = d'^V/dz^ = at the centre of the trap, and 
the electric field Ez{z — 0) along the z axis at the centre of the trap was calculated. The 
voltage on the displaced electrode was then adjusted to null this axial field, and the resulting 
potential shape V{x, y, z) examined. This was repeated for displacements of the electrode in 
the X and y directions. The results are tabulated in table|Hl The voltages are given such that 
for p = 10 /xm and ii^max — 100 V/xm~^, the axial centre of mass frequency is cji = 27r x 22 
MHz, and the radial centre of mass secular frequencies are = 27r x 41 MHz and cjj^ = 27r x 52 
MHz, for calcium ions. 

In this study we did not introduce any further electrodes, whereas in practice further 
electrodes would be present if the structure is part of a larger array, and such further electrodes 
could be used to cancel stray electric fields along z such as the one obtained here. However 
by using the given set of electrodes we obtain an estimate of the precision required for the 
d.c. voltages, and an upper limit on the size of unwanted d.c. quadrupole terms which are 
introduced. 

The values of Zc in table El show that when the voltages are adjusted to cancel the 
quadrupole, then the stray electric field (owing to manufacturing imprecision) is large enough 
to cause a displacement of order p of the ions, which would cause a problem. Comparing the 
values of V2 in the 3rd and last columns of the table, however, it is seen that only a ^ 5% 
adjustment in V2 is needed to cancel this field. Since we are dealing with a 5% misplacement 
of the electrode, we infer that the required values of the electrode voltages are not especially 
sensitive to electrode misplacements. Equally, the figures show that the behaviour is sensitive 
to voltage inaccuracies, because without this same 5% voltage adjustment there would exist 
a severe problem. 

To characterise the latter, argue as follows. In the presence of a field Eq the centre of 
mass of the ions is displaced by approximately z^. — qEo/mujf, and at the octupole condition 
this evaluates to Zc — Eq/ScPP (using eq. Hl^fl). This is a problem when the displacement is 
such that the potential hill introduced when a becomes negative does not appear in between 
the ion pair, so fails to separate the ions. Therefore the displacement is unacceptable if it 
exceeds approximately d/2, so we require 

Eo < Sp. (38) 

In the example under discussion, /3 ~ 2.5 x 10'^^ Vm~^ and d ~ 1 /xm, so the requirement is 
£^0 < 2500 V/m. It is in practice easy to cancel fields of this magnitude, but in view of the 
fact that patch potentials and manufacturing imprecision can result in fields much larger than 
this, it may be necessary to characterise the electrode structure in situ in a vacuum system 
by transporting ions through it. 

Next, consider the d.c. quadrupole terms. When the electrode structure lacks symmetry 
owing to its manufacturing imprecision, we can no longer assume mixed derivatives such as 
d^V / dxdy, d^V / dydx are zero, so we expect to get a d.c. quadrupole in the xy plane when 
the electrodes are adjusted to satisfy the constraints dV/dz = 0, d^V/dz^ = 0. In order to 
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estimate the typical size of d.c. quadrupole caused by manfacturing imprecision, we left the 
electrode voltages at the values indicated in table El without attempting any further nulling 
of quadratic terms, and examined the quadratic terms in the potential near the origin. These 
terms were of order 10^ Vm~^. Under the conditions of the test, the radial confinement 
from the a.c. quadrupole is such as to achieve radial secular frequency around 45 MHz. We 
find that the influence of the d.c. quadrupole is small: it merely changes the radial secular 
frequencies to - (45^ ± lO^)^/^ ~ 45 ± 1 MHz. 

7 Discussion 
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Fig. 19. Values of 7 and /i for the structures considered in the text, (a); high aspect ratio, g > 1; 
*: 1.1-1.3, o: 2.1-2.6, ®- 3.1, O: 4.1, 4.2, 4.7, 4.14. (b): low aspect ratio, g < 0.75; o: 3-layer, o: 
2-layer, +: 'planar' 5.1,5.2 {g ~ 0.65). The filled symbols represent results for g ~ 0.4, the open 
symbols represent results for g ~ 0.2. The dotted lines give the values which were obtained for 
the two 'standard' cases described in section |S] for comparison. 



A summary of the values of 7 and fi obtained for the various structures we have considered 
is given in figure El We find that, with the exception of results 1.1-1.3, 7 is always smaller 
than 11. This may be because we have not allowed the r.f. electrode to be divided. It presents 
a flat potential surface along the z direction which tends to suppress the variation in potential 
which is needed for large 7. It may be possible to increase 7 without introducing multiple r.f. 
electrodes, by tailoring the shape of the edge of the r.f. electrodes to follow an equipotential 
produced by the d.c. octupole. 

One might expect a competition between fi and 7 to be observed in the results, but there 
is little indication of this. The best structures give high 7 and high /i simultaneously. 

Many of the optimizations studied, by adjusting the widths and placements of electrodes, 
gave useful improvements in 7 and/or fi. This implies that a careful optimization of any 
chosen design is worthwhile. 

Of the three-layer electrode structures we examined, those of the general form shown in 
figure H performed best. The highest 7 value (7 ~ 0.03, result 2.2), was obtained with /i close 
to the highest value, using a 3-layer structure of aspect ratio <? — 3. 

At lower aspect ratio, g < 0.5, the two-layer structures performed better than three-layer 
ones of the same g (figure 119b ) . The larger value for 7 may be understood from the fact 
that when g is small, 7 varies rapidly with w because thicker electrodes can extend their 
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influence on V{x, y, z) over larger distances. For a given total depth w of the whole structure, 
if there are m layers and m — 1 gaps all of the same thickness, then the thickness of each 
layer is w/(2m — 1). Therefore a two-layer structure of given g may be expected to give a 
similar value for 7 as a three-layer structure with g larger by a factor 5/3. This may be seen 
by comparing the values for a three-layer structure a.t g = 0.4 in figure El with those for a 
two-layer structure at 17 = 0.24, e.g. result 4.11 (and c.f. figureEb)- The two-layer structures 
also give larger values for fi at given g. This is in part for a similar reason, and also because 
the electrodes are more conveniently placed to realise a radial quadrupole. The two-layer 
structures of large aspect ratio satisfy the condition (|21l) at qr = 0.3 which means that, owing 
to their good radial confinement, it is not necessary to cancel d.c. quadrupole terms in the 
potential, but doing so remains useful in order to reduce the required r.f. voltages. 

At g — 0.24 the highest 7 value (7 ~ 0.005, result 4.11) was smaller than the best overall 
by a factor 6, and fi was reduced by an order of magnitude. 

The planar or near-planar arrangements, cases 5.1 and 5.2, gave a further order of mag- 
nitude reduction in 7 compared to the two-layer low-aspect-ratio structure. Owing to this 
small 7 value, these structures satisfy the condition so the cancellation of d.c. quadrupole 
terms is not necessary. 

8 Conclusions 

We have studied the general problem of separating or bringing together pairs of trapped 
ions, while keeping the frequencies of the normal modes of oscillation of the ions as high as 
possible. A simple picture of this is to say the central electrode introduces a potential hill 
which pushes the ions apart. However we have found that the distance between the ions can 
be small compared to the distance scale of the electrodes, even when the potential well is 
about to divide into two local minima. This means that the potential shape in the region 
between the ions is really a joint property of the whole set of electrodes. There are various 
contributions to the potential that have a quadratic dependence on position near the origin, 
and when the ions are being split or recombined these are balanced, so as to leave a quartic 
dependence. 

The problem under consideration reduces to two main issues: radial confinement, and 
axial confinement at the point where it is weakest, which is to say when the axial potential 
is quartic. We found that axial confinement is the harder to achieve in large structures, and 
radial confinement in small ones. The size scale pc when the radial confinement becomes the 
limiting problem is of order microns. It is set largely by a length scale Lq (equation (|25|l ) 
which is determined purely by materials considerations, i.e. by the maximum electric field 

1 /2 

allowed on electrode surfaces. Since Lq varies as Eninx it will not change greatly for different 
materials or fabrication methods. Ion trap micro-fabrication methods being attempted at 
present are in the regime p ^ pc ? 

We have studied a variety of electrode structures and characterised them mainly through 
the geometrical factors fi and 7 (equations (|20|l and (|19|l 'l which characterise the strength of 
r.f. quadrupole and d.c. octupole moments. Many of the structures exhibit z-dependence 

"Methods have been suggested to approach the nanometre scale, such as the use of carbon nanotubes as 
electrodes. The use of very tight traps would require novel non-optical methods to achieve quantum logic 
gates, because the Lamb-Dickc parameter would be very small. 
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in the oscillating part of the potential, and therefore micro-motion along the z axis. This 
motion is found to have a small Matthieu (/-parameter and is not expected to be a problem? 
We assumed reflection symmetries in the structures, and for one example case we examined 
the effects of manufacturing imprecision which broke the symmetry. This did not suggest any 
major problem with the design, in that the unwanted electric field owing to manufacturing 
error could be cancelled by small adjustments to the electrode voltages, and the unwanted 
quadrupole terms were small. 

We found that for large aspect ratio, 3-layer designs gave the highest values of 7 and 
simultaneously almost the highest /i. At low aspect ratio, 2-layer designs performed best at a 
given aspect ratio. Current efforts to fabricate trap arrays are focussed on distances of order 
p ~ 100 ^m and this suggests a low aspect ratio g ~ 0.2 will be necessary if lithographic 
methods are used. In the future, however, there may be interest in p ~ 10 /im and then g ~ 2 
will be available. 

It is not surprising that the 3-layer designs offer tighter traps than the planar designs, but 
it is noteworthy that when considering the problem of separating ions, this is especially true. 
The factor ~ 150 decrease in 7 (comparing case 2.2 in table|21with case 5.1 in tabled means 
that the planar design would need to be fabricated at a distance scale a factor 150^/^ ~ 5 
times smaller in order to obtain the same motional frequency wi (equation (|22|l '). If we assume 
heating rates scaling as the heating rate in the latter case would be increased compared 
to the former by a factor ~ 800. If instead one compares structures having the same number 
of phonons of heating per ion separation time (c.f. section [3.1(1 . then the planar design leads 
to a motional frequency a factor 150°'^^^ ~ 7 lower. 

In view of the greater ease of fabrication of planar structures, a sevenfold reduction of 
motional frequencies might be regarded as acceptable. However, for a quantum computer in 
which moving ions around and separating and combining strings of ions are essential ingredi- 
ents to most logical operations, the separation/combination time may be the main limitation 
on the overall logic speed. To establish a definite design preference between a small planar 
structure and a larger layered structure, a greater understanding of the heating mechanism 
will be needed. 
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Appendix 1: Axial micromotion 

Consider the axial motion of an ion in the potential given in equation H14() : 

V{0, 0, z) = az^ + Pz^ + (pt^z'^ + lirz'^) cos{nt) (39) 

where f2 is the frequency of the voltage applied to the r.f. electrodes, and we added an 
octupole term to the oscillating part in case that is important (it will turn out that it is not). 
For small departures from equilibrium, the equation of motion for ion i of an n-ion chain 
(assuming the ions are confined strongly in the radial direction) is 

m-^ = ~q{2az, + APzf + {2a, Zi + iPrzf) cos(rjt)) + ^ %(z, - zj) (40) 

j 

where Zi is the position of ion i, and kij is the spring constant arising from the Coulomb term 
in the expansion of the ion-ion repulsion (this expansion has only been taken to first order in 
the ion separation). 

When we sum the equations for each of two ions, and the ion-ion repulsion term drops 
out, leaving 

1^ + ^) = -9(2« + cos(m))(zi + Z2) + (4/3 + 4/3, cos{nt))izl + z^). (41) 

Now substitute Zc = {zi + Z2)/2, and zi — Z2 = d + ( where d is the equilibrium separation of 
the ions and C, represents excursions of the ions from their equilibrium separation. Equation 
(|lT|l becomes 

d'^z 

m-^ = -2q{a + a, cos(rJi))zc - q{P + cos{nt)){3{d + ( f + iz^)zc. (42) 

In the limit (, Zc <^ d we neglect the terms of order (d, and in the second term on 
the right hand side. The resulting equation can be transformed to the standard form of the 
Matthieu equation 

d'^z 

-^ + {a,- 2q, cos(2e))z, = (43) 

using the substitutions 

where we assume the d.c. octupole is such that pd^ a, but the a.c. octupole term can be 
neglected since pr ^ otzja^ <C ctz/d^ . The secular frequency associated with the solutions of 
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equation H43|l is given by 




(45) 



— (2a + 3(3(f + 



2 



) 



(46) 



This equation should be compared with 1)10(1 . We can obtain ((46|l from Hl()|l by the substitution 



The same result is obtained by treating the effect of the oscillating part of V{0, 0, z) using a 
static pseuodopotential. 

In the electrode structures discussed in this paper, the oscillating part of the potential 
has a quadrupole form with axial coefficient always less than half the radial coefficients, 
\az\ < {l/2)\ar\- This means that if we choose the radial Mathieu g-parameter Qr < 0.3, then 
\qz\ < 0.15. In the most extreme case we find that to obtain a' = we require a ~ — 0.04|ar|- 
We find that this requires only a small adjustment to the voltages, compared to the values 
required for a = 0. 

Appendix 2: Basic octupole constructions 

We list some combinations of point and line charges that produce octupoles using a min- 
imal number of separate points or lines. 

8. 1 Systems of point charges 

8 point charges q at the corners of a regular cube of side 2a produce an octupole having 
cubic symmetry (PH]) with (3 — 7g/81-\/37reoa^. 

If we search for a system of 8 point charges all lying in a single plane y = 0, the only 
solution is the one produced by the construction %\2\. 1(33(1 . This is shown in figure EHk . In 
the positive quadrant there is a charge q at (x, z) — (a, d) and a charge —f^q at (/a, fd). For 
example, a — d gives f] — (— 13(7/128\/27reoa^)(l — 1//^); the coefficients of x'^, y'^ and are 
then in the ratios 13 : — 12 : 13, and the potential is a combination of axial and 2D octupoles 
in the ratio —24 : 35. When d = then each set of 4 charges on its own produces a 

potential with d'^V/dz'^{Q, 0,0) 0, and the two sets as described produce equal and opposite 
a2l//ax2(0,0,0). This gives (3 = (-14^273 g/8l7reoa5)(l - l/p); the coefficients of a:^ / 
and are in the ratios 26 : —54 : 56. Similar remarks apply with x and z swapped when 
d = V2a. 

8.2 Systems of ring and line charges 

The potential produced by a semi-infinite straight line charge is V ~ (A/27reo) ln(r + z), 
in a coordinate system such that the line charge lies on the negative z axis and ends at the 
origin (r = (x^ + + z^)^/^). 

An axially symmetric octupole ((29(1 can be created by a single pair of identical ring charges 
on a common axis, whose separation is ^/2 times larger than their radius (shown in figure OUh ). 
For charge per unit length A and radius a, /3 = — 14-\/2/3A/81eoa''. This charge distribution 




a' = a + - \qzaz\ 



(47) 
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can also be interpreted as a rotation about the z axis of point charges at the 'special' location 
X — which we identified in the preceeding section. 

Addition of a further pair of rings allows a solution at any values of the ring separations, 
by adjusting the free parameter of the relative charge per unit length. 

A set of 8 identical line charges, ending at the corners of a regular cube of side 2a, and 
extending outwards along the extended major diagonals, produces an octupole having cubic 
symmetry (13011 with /3 = — 56A/25927reoa^, where A is the charge per unit length. 

Next let us consider two sets of 4 line charges lying in two parallel planes (construction 
(E) of section UTT) . To reduce the number of parameters, we consider semi- infinite lines. We 
take the origin at the centre of symmetry, and orient axes so that the two planes are parallel 
to the X — z plane. Apart from an overall scale factor, there are then three parameters which 
describe the layout: the angular coordinates of the end of one of the lines, which may be 
specified for example by spherical polar angles 0, cj), and the angle ip between the line charge 
and the z-direction. This is shown in figure I20l i. Since we have three parameters {9, cj), ip) 
and two constraints (|35|) . one parameter may be chosen arbitrarily and the others adjusted to 
find a solution. For example in the construction (E) we pick ip and then adjust and (f>. We 
will list some example solutions. For convenience, we specify 9, (j) by giving the rectangular 
coordinates {x, y, z) of the end of the line lying in the positive octant, with an arbitrary scale 
factor. Example solutions are {x,y,z,ip) = (1,1,0,0), (1, 1.434, 0.578, 7r/6), (1, 1.799, 1, 7r/4), 
(1, 2.482, 1.731, 7r/3), (0, 1, l,7r/2). 

A set of 4 parallel infinite line charges passing through the corners of a square in a plane 
perpendicular to them ffigure EOb) produces a 2D octupole (|31l) with /? — X/Sneoa'^, where 
the square has side 2a (see figure EHfeV 

A set of 4 infinite line charges passing along the sides of a cube as in figure produces 
an octupole having /3 — — 3A/47reoa'^ , at the centre of the cube, where the cube has side 2a. 

Next, consider sets of semi-infinite line charges all lying in a single plane. Let the plane 
be the x — z plane and assume symmetry under reflections in the x and z axes. First suppose 
there is only a single line charge in the positive quadrant (so 4 line charges in total). The 
geometry has two parameters: the angle 9 between the z axis and a vector from the origin to 
the end of the line, and the angle ip between the line itself and the z-direction. This suggests it 
might be possible to obtain d'^V/dx'^ = d'^V/dz^ = by a good choice of 61 and tp, but the pair 
of simultaneous equations has no solution. The construction (D) (eq. JSSJ) leads to an infinite 
set of solutions with 8 electrodes in two groups of four. The two groups have the same values 
of 9 and tp; the charge per unit length on the outer group is greater than that on the inner 
group by a factor /^. There are further solutions when the two groups have different angles. 
As an example case, shown in figure EHb. 9 = ir/A, ip — tt (line charges parallel to the x axis 
and finishing on the line x = z) gives /3 = -(15A(16 -I- llV2)/167reo(l + V^^a^Xl - l/P). 
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Fig. 20. Examples of octupolc contructions using line charges. 
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Table 3. Results from numerical calculations for the two-layer electrode structures. The value 
of -Emax calculated using an r.f. electrode voltage of 10® V. Results 4.1 to 4.12 use cuboid 
electrodes. Results 4.13 to 4.17 use cylindrical d.c. electrodes in the arrangement shown in figure 
1131 4.13 to 4.15 are for case (b) in figure ImI 4.16 and 4.17 for for case (c) in figure ImI 



Result 


9 




di 


t 


10*/3 


^max 


7 


-^max 






units 




(P) 


(P) 


(P) 


(Va--) 


(Va-^) 


(10-^) 


(lO'^Va-^) 






4.1 


2 


0.16 


0.39 





-710 


21 


-7.0 


1.77 


0.147 


0.154 


4.2 


2 


0.16 


0.39 


0.39 


-728 


23 


-7.4 


1.70 


0.114 


0.124 


4.3 


2 


0.16 


0.39 


0.78 


-872 


25 


-7.5 


1.53 


0.097 


0.109 


4.4 


2 


0.16 


0.78 





-131 


7.8 


-1.30 


1.32 


0.199 


0.199 


4.5 


1.1 


0.19 


0.94 





-98 


9 


-1.3 


2.25 


0.122 


0.158 


4.6 


1.1 


0.19 


0.47 





-1304 


20 


-7.7 


2.42 


0.122 


0.141 


4.7 


0.4 


0.08 


0.40 





-1563 


34 


-4.6 


4.37 


0.041 


0.042 


4.8 


0.4 


0.08 


0.40 


0.5 


-1848 


33 


-5.5 


4.59 


0.028 


0.031 


4.9 


0.24 


0.08 


0.20 





-1094 


233 


-0.47 


12.6 


0.007 


0.011 


4.10 


0.24 


0.08 


0.40 





-1250 


38 


-3.3 


12.7 


0.011 


0.012 


4.11 


0.24 


0.08 


0.40 


0.5 


-1636 


35 


-4.7 


12.8 


0.0099 


0.011 


4.12 


0.24 


0.08 


0.40 


1 


-1649 


34 


-5.0 


12.7 


0.0061 


0.0068 


4.13 


1.28 


0.24 


0.17 





-399 


229 


-1.6 


19.1 


0.097 


0.100 


4.14 


2.56 


0.42 


0.28 


0.8 


12724 


257 


8.5 


16.8 


0.120 


0.112 


4.15 


2.56 


0.42 


0.4 


0.4 


-482 


9.2 


-9.0 


16.2 


0.167 


0.136 


4.16 


2.56 


0.61 


0.42 


1.21 


-970 


10 


-5.4 


36.5 


0.030 


0.041 


4.17 


2.56 


0.61 


0.42 


0.4 


-885 


10 


-5.0 


33.3 


0.048 


0.067 



Table 4. Voltages on the d.c. electrodes required to produce the results presented in tabic 



Result 


Vi 


V2 


Vs 


units 


(V) 


(V) 


(V) 


4.1 


-1.8451 


1.4055 


-3.2924 


4.2 


-1.6997 


2.2963 


-3.1155 


4.3 


-1.8592 


2.4805 


-3.3316 


4.4 


-0.0885 


0.3648 


-1.6934 


4.5 


-0.0290 


0.2095 


-1.6090 


4.6 


-0.7951 


1.4055 


-2.5750 


4.7 


-1.4109 


2.0817 


-3.1413 


4.8 


-1.4960 


2.1987 


-3.0247 


4.9 


-11.582 


10.442 


-6.3685 


4.10 


-1.3316 


1.9730 


-3.0171 


4.11 


-1.2679 


1.8959 


-2.7263 


4.12 


-1.1612 


1.9083 


-2.5576 


4.13 


-22.030 


4.7513 


-4.5768 


4.14 


25.107 


-5.2939 


18.252 


4.15 


-0.3793 


0.1728 


-1.3831 


4.16 


-0.1479 


0.1103 


-0.8449 


4.17 


-0.1427 


0.0917 


-0.7954 



Table 5. The two cases studied which produced almost-octupoles above the plane of the electrodes. 
Result 5.1 uses the single plane electrode geometry shown in figure lTTI Result 5.2 uses the "railway 
track" geometry shown in figure ITsl 



Result 


Va 


Vb 


Vc 


Ve 


10*/3 


^max 


P 


7 




Py 


units 


(V) 


(V) 


(V) 


(mV) 


(V/a*) 


(V/a) 


(a) 


(10-^) 


(10-^) 


(10-^) 


5.1 


-4.257 


4.560 


-3.527 


-2.197 


-0.0155 


44.74 


0.829 


0.198 


10.4 


10.5 


5.2 


1.089 


1.212 


-4.526 


-201 


0.000247 


20.47 


2.5 


0.189 


3.5 


3.9 
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Table 6. Example of the effect of manufacturing imprecision. The first three rows show the 
voltages required for d'^V/dx'^ = d'^V/dy'^ = d'^Vjdz^ = at the centre of the trap. Ez(z = 0) is 
the axial electric field at z = when this condition is fulfilled. Zc is the displacement of the centre 
of mass of an ion pair in a (l/2)mu;^z^ potential well in the presence of electric field Ez- Vj""" is 
the voltage on the displaced electrode which is required to null the electric field. In all cases the 
r.f. electrode voltage amplitude is ~ 240 V, so that uJr — ■ 



direction 


Vi 


V2 


Vz 


Ez{z = 0) 


Zc = qEz/muj^ 






(V) 


(V) 


(V) 


(V/im-^) 


(^tm) 


(V) 


X 


281.7 


97.6 


593.0 


-0.17 


21 


85.8 


y 


281.6 


100.7 


593.0 


-0.058 


7.3 


97.4 


z 


277.3 


100.0 


593.0 


-0.073 


9.2 


95.1 



